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Abstract 

Aims. The Earth Occultation Technique (EOT) has been applied to Fermi's Gamma-ray Burst Monitor (GBM) to perform all-sky 

monitoring for a predetermined catalog of hard X-ray /soft y-ray sources. Imaging with a Differential filter using the Earth Occultation 

Method (IDEOM) has been developed to search for sources not in the catalog, thus completing the catalog and reducing a source of 

systematic error in EOT. 

Metliods. IDEOM is a tomographic imaging method that takes advantage of the orbital precession of the Fermi satellite. Using 

IDEOM, all-sky images have been generated for ~ 4 years of GBM data in the 12-50 keV, 50-100 keV and 100-300 keV energy bands 

in search of sources otherwise unmodeled by the GBM occultation analysis. 

Results. Analysis resulted in the detection of 43 sources in the 12-50 keV energy band, 23 sources in the 50-100 keV energy band, 

and 7 sources in the 100-300 keV energy band. IDEOM analysis has resulted in the addition of 16 sources to the GBM-EOT catalog. 

We also present the first joined averaged spectra for four persistent sources detected by GBM using EOT and by the Large Area 

Telescope (LAT) on Fermi: NGC 1275, 3C 273, Cen A, and the Crab. 

Key words. Keywords should be given 



1. Introduction 

At hard X-ray/soft gamma-ray energies, sources vary on 
timescales from fractions of seconds to years. In order to ob- 
serve transient emission, instruments must be capable of con- 
tinuously monitoring the sky with a wide field of view. At 
present, the current wide-field hard X-ray monitoring instru- 
ments are the Gas Slit Camera (GSC) on the Monitor of All- 
sky X -ray Imaging (MAXI) , which covers the 2 to 30 keV re- 
gion (iMatsuoka et alj l2009l) . the Burst Alert Telescope (BAT) 
on Swift, which covers the 15 to 150 keV band (Gehrels 2004), 
and the Imager on Board the INTEGRAL Satellite (IBIS), which 
covers the 15 keV to 10 MeV range (Win kler et al. 2003) . 
There has not been an all-sky monitor observing low energy 
gamma rays above ~ 200 keV since the Burst And Transient 
Source Experiment (BATSE) on board the Compton Gamma Ray 
Observatory (C GRO) from 1 991-2000, which covered the range 
20 to 1000 keV (iHarmon et al...2004.) . 

In 2008 June the Fermi satellite was launched with two wide- 
field of view instruments: the Large Area Telescope (LAT) 
and the Gamma-ray Burst Monitor (GBM), spanning the ~ 
20 MeV to 300 GeV ( Atwood et al. 2009) and 8 keV to 40 MeV 
dMeeean et aT]l20()9l) energy bands, respectively. GBM is cur- 
rently the only instrument performing all-sky observations above 
~ 200 keV and below the 20 MeV LAT threshold. The Earth 
Occultation Technique (EOT) has been applied to the GBM 
detectors in order to mon itor fluxes of known point sources 
dWilson-Hodge et all l2012h . Over ~ 4 years of Earth occul- 
tation analysis, the GBM Earth Occultation team has gener- 



ated daily updated light curves for the ~ 200 sources cur- 
rently monitored (available at http://heastro.phys.lsu.edu/gbm/) 
and published a catalog of high energy (> 100 keV) sources 
(Case et al. '2 0111). and a ful l 3-year catalog of all GBM sources 
(Wilson-Hodg e et all l2012h . A particularly surprising result 
from GBM's continuous all-sky monitoring was the discovery 
of the variability of the ha rd X-ray flux from the Crab Nebula 
(IWilson-Hodee et aI.ll20Tlh . 

The standard EOT approach requires an input catalog of source 
positions to monitor sources. An additional method is necessary 
for finding sources not in the input catalog. In order to locate 
these sources with the uncollimated GBM de tectors, a tomo - 
graphic imaging method has been developed ( Rodiet £01201 ih . 
This technique employs projections of the the Earth's limb on the 
sky over the course of Fermi's orbital precession period (~ 53 
days) to generate all-sky images. Imaging with a Differential 
filter using the Earth Occultation Method (IDEOM) allows for 
searching for undiscovered sources and known sources absent 
from the GBM input catalog. Having a source catalog that is as 
complete as possible helps reduce systematic errors with EOT 
and may make it possible to investigate the discrepancies seen at 
high energies between the Marshall Space Flight Center (MSEC) 
and the Jet Propulsion Laboratory (JPL ) analyses of BATSE data 
(iHarmon et ; 



Propulsion Laboratory (JPL 
: al.l l2004l:lLing et al.ll2000l) . 



In this paper we describe the IDEOM imaging method and its 
application to GBM, and present all-sky images for the first 
~ 4 years of data in the 12-50 keV, 50-100 keV, and 100-300 
keV energy bands. In Section 2 we give a brief review of the 
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GBM instrument and the standard EOT analysis; Section 3 de- Table 1. Example of systematic error due to uncatalogued 
tails IDEOM and its application to GBM; Section 4 presents the sources 
imaging results. Section 5 discusses the four persistent sources 
detected by GBM and LAT, and in Section 6 we discuss the con- 
clusions and future work. 



2. Fermi's Gamma-ray Burst Monitor 

Fermi was launched in 2008 June to an altitude of 565 km and an 
orbital inclination of 25.6°. At this altitude the Earth's diameter 
is ~ 140° and thus ~ 30% of the sky is blocked by the Earth 
at any given time. Over 85% of the sky is occulted by the Earth 
over the course of a single orbit, and the entire sky is occulted in 
~ 26 days (half of the Fermi orbital precession period), though 
the exposure is not uniform. 

GBM consists of 1 2 Nal detectors and two BGO detectors 
(iMeegan et alJl2009l) . The Nal detectors (covering the range 8 
keV to 1 MeV) are 12.7 cm in diameter and 1.27 cm thick and 
are arranged in four groups of three detectors, located on the cor- 
ners of the spacecraft. Six detectors are oriented perpendicular to 
the z-axis of the spacecraft. Four are inclined +45° from the +z- 
direction, and two are inclined +20° from the +z-direction. The 
BGO detectors (covering energies ~ 200 keV to ~ 40 MeV) 
are 12.7 cm in diameter and 12.7 cm thick, and are located on 
opposite sides of the spacecraft. Since 2008, GBM has had two 
continuous data types, CTIME and CSPEC. CTIME has nom- 
inal 0.256-second resolution and 8-channel spectral resolution. 
CSPEC has nominal 4.096-second resolution and 128-channel 
spectral resolution. 

As shown with BATSE on CGRO, monitoring known hard X- 
ray/soft gamma-ray sources is possible with simple non-im aging 
detectors with EOT (iLing et al.ll2000l; iHarmon et al.ll2002 1 As a 
source is occulted by the Earth, the detector count rate decreases 
like a step function. A similar feature, but with opposite polar- 
ity, occurs in the detector count rate as the source comes out of 
occultation. Occultation times can be predicted using the source 
coordinates and the spacecraft positions. The occultation time 
is defined to be the time at which the transmission through the 
atmosphere is 50% at 100 keV. Given a predetermined catalog 
of sources, the flux for a source of interest is calculated by se- 
lecting a 4-minute window of count-rate data centered on the 
occultation time and fit to a quadratic background plus a source 
transmission model for each source that occults during the win- 
dow. A source transmission model is generated that consists of 
an assumed source spectrum combined with an energy depen- 
dent atmospheric transmission model, which is then convolved 
with the detector response. For each energy channel for each 
detector, a scale factor is derived from the fit of the source trans- 
mission models to the data and averaged over all of the detectors 
viewing the source. The weighted average scale factor for the 
source of interest is multiplied by the spectral model flux to de- 
termine the final flux. The EOT applied to GBM is described in 
greater detail in Wilson-Hodge et al. (2012:) . 



3. Earth Occultation Imaging 

EOT can monitor only those sources included in the input cat- 
alog. For the CGRO/BAlTSE mission, two tomographic imag- 
ing methods were developed. E arth occultation transform imag- 
ing dZhang et all 1 19931 Il995h and the Likelihood Imaging 



Method for BATSE Occultation data (LIMBO) dShaw et alj 
2004). Earth occultation transform imaging generated images 
with the BATSE data using an inverse Radon transform and 
a maximum entropy method (MEM) deconvolution. With this 
technique, images were created for time periods spanning days 
to weeks and energies from 20 to 300 keV. The size of these 
images ranged from 5° x 5° to 40° x 40°. Regular analysis con- 
sisted of imaging the galactic plane with 27 30° x 30° fields of 
view to search for sources. This method was a ble to discover a 
number of transient sources (Zh ang et al.lll995l) . but has charac- 
teristics that make systematically searching for sources difficult. 
First, the effective image size is limited because of distortions 
due to approximating the Faith's limb as a straight line instead 
of as an arc, resulting in inaccurate source locations. A second 
difficulty is caused by the projections along the limb of bright 
sources, which can extend for tens of degrees. Earth occultation 
transform imaging has no way to account for these projections 
which bias the image background. Using MEM, which is an iter- 
ative, non-linear process, that results in relative source intensities 
that include spurious contributions of other sources in the field 
of view, which makes comparing different images difficult. 

LIMBO generated all-sky images using a likelihood statistic 
for a giid of predefined source positions on the sky. In this 
method the data were background subtracted using a satellite 
mass model and then passed through a differential filter before 
calculating the Maximum Likelihood Ratio for that location. 
Sources were then found using a CLEAN deconvolution algo- 
rithm. All-sky images were generated for 489 days of data for a 
single broad energy band covering 25-160 keV. While LIMBO 
exhibited improvements over the Earth occultation transform 
technique, its resolution was limited due to the 2° separation of 
its sky grid points. Also by using only one broad energy bin, 
spectral information about sources was lost. 

The cuiTent IDEOM approach makes it possible to locate un- 
modeled sources by imaging the entire sky. Compiling as com- 
plete a catalog as possible will reduce a source of systematic 
errors in the EOT analysis which may arise by attributing flux 
from an unmodeled source to a source actually in the catalog. 
This problem can been seen for example with IE 1740-29 and 
GX 1+4 (separated by ~ 5.5°), which occasionally goes into 
outburst. If GX 1 +4 is absent from the catalog, its transient flux 
is attributed to nearby sources (e.g. IE 1740-29). The first two 
columns of Table [T] show the significance of IE 1740-29 and 
GX 1+4 from MJD 55708 to 55774. During this time GX 1+4 is 
flaring. The third column shows the significance for IE 1740-29 
during the same period when GX 1 +4 is omitted from the cat- 
alog. The table shows an increase of ~ 3cr in each energy band 
for IE 1740-29 when GX 1+4 is not taken into account. Thus 
having a catalog that is as complete as possible is necessary for 
accurate flux measurements. 

Another motivation is to investigate a potential cause of the dis- 
crepancies between the MSEC and JPL analyses of the BATSE 
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occultation data at high energies (100-300 keV). Results from 
the JPL method (iLing et alJl2000l) have shown larger high en- 
ergy fluxes for som e sources than corre sponding results from 
the MSFC analysis dHarmon et alJl2004l) . These apparent high 
energy tails appear correlated with the phase in the 52-day pre- 
cession period of the satellite, indicating a systematic effect pos- 
sibly due to unmodeled sources. 



3.1. IDEOM with GBM 

For a single occultation, the change in count rate corresponds 
to a projection of the brightness of a source along the Earth's 
limb at the time of the occultation. A fixed location on the sky 
is occulted at different points along the limb on successive or- 
bits. Each point along the limb corresponds to a different limb 
projection angle. Thus during the orbital precession of the satel- 
lite, a range of projection angles are sampled. These projections 
can then be combined to generate an image of the sky with an- 
gular resolution dependent on the duration of the occultation 
step. In the energy range corresponding to the GBM Nal de- 
tectors (8 keV to 1 MeV), a single step has a finite duration 
of ~ 8s/ cos/3 due to absorption by the Earth's atmosphere, 
where /3 is the elevation ang le of the occult e d source with respect 
to th e sateflite's orbit (iCase et al.l 120111: IWilson-Hodge et al.1 
I2OI2I) . The best possible position resolution can be expressed 



as A0 = 360° X 



8 .s / cos/j 



^ where P is the satellite's orbital period 

dShaw et al.ll2004t) . FennVs orbital period of ~ 96 minutes cor- 
responds to a minimum A6 ~ 0.5°. A6 increases as /3 increases 
until y6 > 70°, at which point sources are no longer occulted. 

To generate all-sky images using occultation imaging, an input 
catalog is created of ~ 600, 000 virtual sources uniformly cover- 
ing the entire sky with 0.25° spacing. Energy dependent occul- 
tation times are predicted for all the virtual sources so that the 
occultation time occurs at 50% transmission in a given energy 
band. A 4-minute window of GBM CTIME data is selected for 
the detectors that view each occultation at less than 75° from the 
detector normal. A first order detector response correction is ap- 
plied to account for detector efficiency as a function of the pho- 
ton energy and the angle between the source and detector nor- 
mal. Combining both rise and set occultation step types provides 
observations at more projection angles. In order to combine rise 
steps and set steps, windows with rise steps are converted to set 
steps by rotating the window about the occultation time at the 
center of the window. The weighted average of all windows for 
each individual virtual source is calculated over an orbital pre- 
cession period. Poisson errors are assumed for combining oc- 
cultation steps. Error values are then propagated throughout the 
rest of the analysis. This average window is rebinned from the 
standard CTIME 0.256-second resolution to 2.048-second reso- 
lution a nd then pass e d thro ugh a differential filter (See FigHJ as 
used by I Shaw et all (l2004l) . such that the sum of ft, bins to the 
right of the central gap is subtracted from the sum of fi, bins to 
the left with a gap !/„ bins wide between the two sums, and the 
difference is then divided by fh, i.e. 



O: = 






i=i+fc+fb , 



i+f„ 



yj=i-fa 

^H-L-fb 



fb 



(1) 



where rj is the number of counts in bin j and o, is the differ- 
ence between the two sums divided by fi, for a given bin. To 
maximize sensitivity the inner bound should be large enough to 
that 2fa is roughly equal to the occultation duration thus max- 




Figurel. Top: Averaged window for 11 days for the Crab. 
Bottom: The filtered data with the central portion fit to a poly- 
nomial (green) stars, and the outer background portion fit to a 
spline function (red) triangles. (A color version of this figure is 
available in the online journal.) 



imizing the difference between total transmission and total at- 
tenuation. Increasing the outer bound, fh, decreases the angular 
resolution while decreasing ft increases the angular resolution, 
but decreasing ft too much makes the sum susceptible to back- 
ground fluctuations. In practice, we use /„ = 3 and /^ = 8 in 
order to minimize effects from fluctuations and maximize sensi- 
tivity. 

Cuts are made on the raw data to exclude times when the space- 
craft has a high spin rate and when it passes through the South 
Atlantic Anomaly (SAA). High spin-rate effects occur when the 
spacecraft slews, which leads to a rapidly varying background 
during the 4-minute window. When passed through the differ- 
ential filter, a rapidly varying background can mimic an occulta- 
tion step, creating artifacts in the image. The SAA passage effect 
arises from activation in the detectors during passes through the 
SAA and results in large broad spikes in the data, which again 
leads to a rapidly varying background. 

Filtering produces a dip at the occultation time while smoothing 
out background features, as shown in Figure [T]for the Crab. The 
source intensity is calculated from the amplitude of the dip by 
fitting bins within ±{2fj + fb) of the occultation time to a poly- 
nomial, and bins outside this region to a spline function joined 
by a straight line. The amplitude of the virtual source is found 
by taking the difference between the two fits at the occultation 
time. A region of +(2fa + fb) is used to ensure that the dip is 
fully enclosed and thus the difference between the dip and the 
background is optimally sampled. The bottom panel of Figure [1] 
shows the filtered window with the polynomial fit in (green) stars 
and the spline fit in (red) triangles. The outer dashed lines are the 
boundaries for the polynomial fit, and the middle dashed line de- 
notes the occultation time. 
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3.2. Systematic Effects 

Two systematic effects are explicitly taken into account in the 
processing: 1) source confusion with bright sources and 2) win- 
dows in which the rise and set steps of a bright source dominate 
the spline background fit. The issue of source confusion arises 
from the ambiguity of where along the Earth's limb the measured 
flux originates. Since the GBM detectors have no direct imaging 
capability, the flux for an occultation time can be attributed to 
any location on the sky along the Earth's limb. Consequently, 
bright sources (e.g. Crab, Cyg X-1, Sco X-1) produce a charac- 
teristic "X" pattern that can extend for tens of degrees on the sky. 
These "arms" are the result of projections remaining in a simi- 
lar orientation for a relatively long time as the source reaches its 
maximum and minimum /? angles during the Fermi precession 
period, and may ovei-whelm any flux that comes from possible 
faint sources lying on or near one of these arms. To mitigate 
the effects of source confusion, an algorithm has been developed 
to ignore occultation windows for virtual sources when bright 
sources occult close in time (< 1 1 seconds) but far away on the 
sky ( > 10°). 

The other systematic effect involves virtual sources at roughly 
the same declination as a bright source and a few degrees away 
in right ascension. The virtual source is close enough on the sky 
to the bright source so that both rise and set steps for the bright 
source occur inside the 4-minute window, but the two sources 
are far enough away such that when the window is filtered, the 
rise steps and the set steps for the bright source remain sepa- 
rated. Consequently, the steps from the bright source are treated 
as the background level. Since the background level may then 
be larger than the source at the occultation time, the calculated 
amplitude of the virtual source appears negative. To reduce this 
effect, the polynomial fit boundaries are shifted towards the oc- 
cultation time (i.e. the center of the window) to reduce the effect 
of the bright source on the background fit. 



Table 2. Fermi Orbital Precession Periods 




40 60 

of 2-sec time bins 



Filtered Averaged Window for Crab for Aug 12, 2008 to Got 4, 2008 




of 2-sec time bins 



Figure 2. Top: Averaged window for 1 day for a virtual source 
near the Crab. Crab steps at bins are ~ 35 and ~ 75. Bottom: The 
filtered data with the central portion fit to a polynomial (green), 
and the outer background portion fit to a spline function (red). 
Large dips from Crab steps cause calculated background level 
to be significantly larger in magnitude than center of window, 
resulting in negative amplitude for virtual source. 
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Fig.|2]demonstrates this problem by using a virtual source at a = 
80.75°, (5 = 22.0°, close to the Crab (a = 83.63°, 6 = 22.02°). 
The top panel in Figure 2 shows the rise and set occultation steps 
for the Crab, for one precession period, roughly at bin 35 and bin 
75. The bottom panel shows how the filtered Crab steps misrep- 
resent the background and thus the calculated amplitude of the 
virtual source. 



4. Results 

IDEOM has been applied to the first ~ 4 years of GBM CTIME 
data, covering 2008 August 12 to 2012 September 14, which 
corresponds to 28 Fermi orbital precession periods. An all-sky 
image has been generated for each precession period for three 
broad energy bands. Table |2] shows the beginning and ending 
dates for each precession period. Single precession period im- 
ages for an energy band were combined by summing the inten- 
sities for a virtual source weighted by the number of occulta- 
tion steps for the point during a precession period. Figs. |3]|5] 
show the combined significance maps for the first ~ 4 years of 
the mission for the 12-50 keV, 50-100 keV, and 100-300 keV 
energy bands, respectively. Contours have been over plotted at 
3.5cr, 5cr, lOcr, and 20cr. The source names plotted in Figs. 3-6 
correspond to sources detected > 5cr during the 28 precession 
periods. Precession period 22 was excluded from the maps of 
the lowest energy band due to high significance fluctuations in 
the background. Also, time periods when the high mass X-ray 
binary A0535-H262 was in an extremely bright outburst, during 
which it flickered, resulting in wave-like features contaminat- 
ing the images were also removed. These times occurred during 
precession periods 9, 10, and 18. Early IDEOM analysis con- 
sisted of cross-correlating features in images with Swift/BAil, 
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INTEGRAL/SPl, and FermifLA^ catalogs, and then analyzing 
that source with EOT to determine fluxes and significances for 
that source. Roughly 15 sources were added to the GBM-EOT 
catalog this way. 

To determine positions for potential sources in the images, the 
sky was searched for points with intensities greater than 3 .5cr 
above background. 5° x 5° regions centered on each point were 
fit to a two-dimensional Gaussian. The maximum value of the 
fit in standard deviations above the background was taken as the 
source significance, and the centroid was taken to be the source 
position. Using the parameters of the fit, the integrated signif- 
icance within a rectangle ±2<Xa,6 of the source position is cal- 
culated, where 0-^ and os are the calculated width of the fitted 
Gaussian. 

Candidate sources were selected from the list of potential 
sources by requiring that the product of the peak significance 
and the integrated significance divided by the statistical error for 
the virtual source be > 10. This statistic is to ensure that vir- 
tual source selection was not biased by individual "hot" pixels. 
Candidate sources for an image are used in a temporary input 
catalog and analyzed using the standard EOT to determine fluxes 
and significances. The cumulative images for each energy band 
are searched to look for candidate sources. Tables |3]- |5]contain 
the list of sources identified by IDEOM and verified by EOT. In 
Fig. [3] the sources with blue triangles denote sources added to 
the input catalog based on imaging analysis. 



4.1. Comparison with GBM EOT Results 

The sources detected by IDEOM in the 12-50 keV band are listed 
in Table |3] In this energy range, the EOT detects 67 sources 
while IDEOM finds 43, two of which are not detected by the 
EOT. Of the 26 sources not found by IDEOM, 1 1 were close to a 
brighter source (within a few degrees) and thus could not be sep- 
arated, 12 did not show a feature above 3.5cr in the image, two 
did not meet the detection threshold of > 10 for the test statis- 
tic, and one (A0535-I-262) lacked image data during outburst 
and was not expected to be detected. The EOT detected 17 tran- 
sient sources in this energy range; 9 of which were detected by 
IDEOM. A sensitivity Hmit for IDEOM in the 12-50 keV band 
can be estimated based on the EOT fluxes for the sources where 
no feature was seen in the image. There were 10 sources with 
an average flux < 25 mCrab detected by the EOT. Only three of 
those were found with IDEOM, suggesting a sensitivity limit at 
~ 25 mCrab for isolated sources. For the 18 sources with aver- 
age fluxes between 25 mCrab and 50 mCrab, six sources were 
missed by IDEOM. These sources are all located in low sensi- 
tivity regions. Two (4U 1702-429 and 4U 1705-472) lie close to 
a Sco X-1 arm thus decreasing the number of usable occultation 
steps. The other four sources are farther away from the Sco X- 1 
arm, but they are at lower declinations so are not occulted during 
part of a precession period and also lose some occultation steps 
due to interference from Sco X-1. Thus for sources that are not 
isolated, the sensitivity limit is closer to 50 mCrab For the two 
GBM sources found by IDEOM, 3C 273 and GRS 1724-308 
have an EOT significances of 4.8cr and 4.7cr, respectively. The 
corresponding IDEOM significances are 5. Oct and 13.5cr. The 
EOT and IDEOM fluxes for 3C 273 are consistent within the 
eiTors. For GRS 1724-308, the IDEOM flux is 3.5 times larger 
than the EOT flux, likely due to contributions from source in- 
terference in the crowded galactic center region. Using the posi- 
tion error between the measured source position and the known 



source position, the positional accuracy in the 12-50 keV band 
is Ae ~ 0.40° 

Table |4] shows the sources detected by IDEOM in the 50-100 
keV band. In this energy range GBM-EOT detects 25 sources 
while IDEOM detects 23 sources, three of which are in the GBM 
catalog but not detected. There were 6 transients detected by the 
EOT with 5 of the found by IDEOM. Three of the five sources 
not detected by IDEOM have no feature in the image (NGC 
2110, Sco X-1, and GX 354-0). The other two sources were close 
to a detected source and unable to be separated by IDEOM (OAO 
1657-415 and lA 1742-294). The average flux of the three un- 
detected sources without an image feature is less than 22 mCrab, 
suggesting that the sensitivity in the 50-100 keV energy band is 
roughly 22 mCrab. The three GBM sources detected by IDEOM, 
but not EOT were NGC 4388, GRS 1758-258, and 4U 1812- 
12. The EOT sensitivity is low for sources close together as 
the number of usable occultation steps for each source is re- 
duced. IDEOM does not exclude occultation steps for sources 
close together on the sky when generating an image allowing 
GRS 1758-258 (~ 0.75° away from GX 5-1) and 4U1812-12 
(~ 1.9° away from GX 17-h2) to be detected by IDEOM and 
not by the EOT. INTEGRAL observations have shown spectra of 
GRS 1758-258 and 4U 1812-12 with signifi cant flux out to 100 
keV dPottschmidt et al.ll2008t 'Tarana et al.'l2006 ^ that the other 
two sources do not posse ss (Paizisetal. 2005: iMainardi et aP 
I2OIO ': 'Miglia ri et aDl2007h . which is in agreement with results 
from IDEOM. The third source, NGC 4388 is detected at 5.4cr 
by IDEOM and 4.7cr by the EOT and with flux measurements 
consistent within the eiTors between the two methods. The cal- 
culated position accuracy in this energy band is A0 ~ 0.47°. 

Table|5]shows the sources found by IDEOM in the 100-300 keV 
band. In this energy range EOT detected six sources. No tran- 
sient sources were detected by EOT or IDEOM in this energy 
band. IDEOM found all six of these sources and also detected 
GRS 1758-258. The lowest average flux is ~ 54 mCrab for 
GRS 1758-258 at 9.0cr. The EOT significance for this source 
is 1.6cr with a flux of 58 mCrab. The difference in signifi- 
cance is due to source crowding and thus fewer occultation mea- 
surements. There are two sources detected by IDEOM between 
3.5cr and 5cr. The sources are 3C 273 with an average flux of 
~ 20 mCrab and a significance of 3.9cr and XTE J 1752-223 
with an average flux of ~ 32 mCrab and a significance of 4.6cr. 
Using XTE J1752-223, the 4-year sensitivity can be estimated 
to be between 32 and 54 mCrab. The lack of unmodeled sources 
at these energies suggests that the catalog is complete down to 
roughly 20 mCrab in the 100-300 keV band and that the dis- 
crepancies between the EOT and EBOP analyses of BATSE oc- 
cultation data are likely due to inaccuracies in the EBOP back- 
ground model. The calculated position accuracy in the 100-300 
keV band is A9 ~ 0.25°. 

The sensitivity difference between the EOT and IDEOM is likely 
due to several factors. Fkst the EOT uses the exact source posi- 
tion and thus is able to optimally sample the occultation step 
(and flux) from a source while the IDEOM virtual source is not 
necessarily close to the source location and not well sample the 
each occultation step. Also, the interference time is larger for 
IDEOM (11 sec) compared to the EOT (8 sec) which can re- 
sult in fewer usable occultation steps for some source locations 
on the sky (e.g. near Sco X-1 arms). Finally, in order to reduce 
potential problems after SAA passes, the data less than 400 sec 
after an SAA pass have been ignored likely losing valid occulta- 
tion steps. 
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Table 3. Sources Detected by IDEOM in the 50-100 keV Band 



Measured Position 


Source Name 


Known Position 


Position Error 


(Degrees) 




(Degrees) 


(Degrees) 


(19.37, -73.50) 


SMC X-1 


(19.27, -73.43) 


0.07 


(58.92, 31.00) 


XPER 


(58.85, 31.05) 


0.08 


(83.59, 22.00) 


CRAB 


(83.63, 22.01) 


0.05 


(83.59, -66.00) 


LMC X-4 


(83.20, -66.36) 


0.40 


(94.16,9.00) 


4U 0614+09 


(94.28,9.13) 


0.18 


(135.63, -40.75) 


VELA X-1 


(135.52, -40.55) 


0.21 


(170.59,-60.50) 


CEN X-3 


(170.31,-60.62) 


0.18 


(176.79,-62.00) 


IE 1145.1-6141 


(176.86,-61.95) 


0.06 


(182.40, 39.25) 


NGC4151 


(182.65, 39.41) 


0.25 


(186.68,-63.00) 


GX 301-2 


(186.65, -62.77) 


0.23 


(187.34, 1.75) 


3C273 


(187.27,2.05) 


0.31 


(195.43,-61.50) 


GX 304-1 


(195.32, -61.60) 


0.11 


(201.27,-43.50) 


CEN A 


(201.36,-43.01) 


0.49 


(244.99,-16.00) 


SCO X-1 


(244.98, -15.64) 


0.36 


(248.25, -67.50) 


4U 1626-67 


(248.07, -67.46) 


0.08 


(254.57, 35.50) 


HER X-1 


(254.45, 35.34) 


0.18 


(255.78, -40.50) 


OAO 1657-415 


(255.19,-41.67) 


1.25 


(256.06, -37.00) 


GX 349+2 


(256.45, -36.41) 


0.66 


(256.39, -39.00) 


4U 1700-377 


(255.98, -37.84) 


1.20 


(261.81,-30.75) 


GRS 1724-308 


(261.89,-30.80) 


0.09 


(262.88, -25.00) 


GXl+4 


(263.00, -24.74) 


0.28 


(263.02, -17.75) 


GX9+9 


(262.93,-16.96) 


0.79 


(263.36, -32.75) 


GX 354-0 


(263.35, -33.39) 


0.64 


(264.41, -43.75) 


X1735-444 


(264.74, -44.45) 


0.74 


(266.55, -28.50) 


MAXIJ1745-288 


(266.46, -28.82) 


0.33 


(266.74, -30.00) 


IE 1740-29 


(265.97, -29.74) 


0.72 


(266.82,-31.50) 


IGRJ17464-3213 


(266.57, -32.24) 


0.77 


(268.71, -2.25) 


SWIFT J1753.5-0127 


(268.36,-1.45) 


0.87 


(270.33, -25.00) 


GX5-1 


(270.27, -25.08) 


0.10 


(270.57, -20.00) 


GX9+1 


(270.38, -20.53) 


0.56 


(273.48, -12.25) 


4U 1812-12 


(273.80,-12.08) 


0.35 


(273.82, -17.25) 


GX 13+1 


(273.63,-17.15) 


0.20 


(273.93, -14.25) 


GX 17+2 


(274.00,-14.03) 


0.23 


(275.93, -30.50) 


HI 820-303 


(275.92, -30.36) 


0.13 


(276.83, -36.25) 


3 A 1822-371 


(276.44,-37.10) 


0.91 


(277.29, -23.00) 


GS 1826-238 


(277.36, -23.79) 


0.80 


(280.02, 4.25) 


SERX-1 


(279.99, 5.03) 


0.79 


(285.06, -25.75) 


HT 1900.1-2455 


(285.03, -24.92) 


0.83 


(288.81, 11.00) 


GRS 1915+105 


(288.82, 10.97) 


0.04 


(299.70, 35.25) 


CYGX-1 


(299.59, 35.20) 


0.10 


(308.07,41.00) 


CYG X-3 


(308.10,40.95) 


0.05 


(308.19,37.50) 


EXO 2030+375 


(308.06, 37.63) 


0.17 


(325.98, 38.25) 


CYG X-2 


(326.17, 38.32) 


0.16 



5. Sources Detected by GBM and LAT 



Three persistent sources are found by IDEOM that are also de- 
tected by the Fenni/LAT. the Crab, 3C 273, and Cen A. A 
fourth source, NGC 1275, is not found by IDEOM, but is de- 
tected by GBM-EOT. Spectra are shown in Fig.|6]for these four 
sources using 4 years of broa d band GBM d ata, LAT 2FGL 
catalog data (Nolan et al."2012')), INTEGRAL (Eckert & Paltani' 
[2009; Jourdain & Roques 2009) data, and data from COMPTEL 
dSchonfelder et al.ll2()00 ). The GBM data span 2008 August 12 
to 2012 August 11, LAT data cover 2008 August 4 to 2010 
August 1, and COMPTEL data are from 1991 May 16 to 1996 
October 15. INTEGRAL observations for NGC 1275 are from 
2003 March to 2004 August while observations of the Crab 
are over the span of ~ 5.5 years from 2003 February to 2008 
September 



5.1. NGC 1275 

NGC 1275 is an elliptical galaxy in the middle of the Perseus 
galaxy cluster. HEAO 1 observations showed emi ssion at ener- 
gies f rom 10-93 keV with no obvious variability (iPrimini et alJ 
119811) . Later observations with CGRO/OSSE failed to detect 
NGC 1275 in the 50-500 keV energy range with upper limits 
significantly be low the fluxes observed during earlier missions 
(IOsakoetalJ1994) . indicating long-term variabihty. INTEGRAL 
has detected the source a t 3-20 keV with JEM- X and at 20- 
120 keV with IBIS/ISGRI (" Eckert & Paltanill2009h . The flux ob- 
served by INTEGRAL was time variable but lower by a factor 
of 3 than the extrapola tion of the £"^ pow e r law component 
report ed by Chandra dSanders et alJ l2005h . lEckert & Paltanil 
(I2009h fit the INTEGRAL data from 2003 March to 2004 August 
with a model of a two-temperature plasma with a central tem- 
perature kT - 3 keV plus an AGN contribution. At higher en- 
ergi es, CGROfEGRET observations resulted in only upper lim- 
its (IReimer et al.ll2003h while the FennifL AT detects NGC 1275 
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Table 4. Sources Detected by IDEOM in the 50-100 keV Band 



Measured Position 
(Degrees) 



Source Name 



Known Position 
(Degrees) 



Position Error 
(Degrees) 



(58.92, 31.00) 


XPER 


(58.85,31.05) 


0.08 


(83.59, 22.00) 


CRAB 


(83.63, 22.01) 


0.05 


(94.42, 8.00) 


4U 0614+09 


(94.28,9.13) 


1.15 


(135.30, -40.75) 


VELAX-1 


(135.52, -40.55) 


0.26 


(182.40, 39.25) 


NGC4151 


(182.65, 39.41) 


0.25 


(185.86, 10.50) 


NGC 4388 


(186.45, 12.66) 


2.24 


(187.61, 2.00) 


3C273 


(187.27,2.05) 


0.34 


(201.48,-43.25) 


CENA 


(201.36,-43.01) 


0.25 


(255.56, -48.75) 


GX 339-4 


(255.70, -48.78) 


0.10 


(255.88, -38.50) 


4U 1700-377 


(255.98, -37.84) 


0.66 


(262.90, -24.75) 


GXl+4 


(263.00, -24.74) 


0.10 


(265.96, -30.25) 


IE 1740-29 


(265.96, -29.74) 


0.51 


(266.07,-31.75) 


IGRJ17464-3213 


(266.57, -32.24) 


0.65 


(266.49, -28.25) 


MAXIJ1745-288 


(266.46, -28.82) 


0.57 


(268.06, -23.00) 


XTE J1752-223 


(268.04, -22.32) 


0.68 


(268.34, -1.50) 


SWIFT J1753.5-0127 


(268.36,-1.45) 


0.06 


(270.36, -26.00) 


GRS 1758-258 


(270.30, -25.73) 


0.27 


(273.77,-13.00) 


4U 1812-12 


(273.80, -12.08) 


0.92 


(274.03, -14.75) 


GX 17+2 


(274.00,-14.03) 


0.71 


(277.49, -24.00) 


GS 1826-238 


(277.36, -23.79) 


0.23 


(288.82, 10.75) 


GRS 1915+105 


(288.82, 10.97) 


0.22 


(299.70, 35.00) 


CYGX-I 


(299.59, 35.20) 


0.22 


(307.89, 40.75) 


CYG X-3 


(308.10,40.95) 


0.26 



Table 5. Sources Detected by IDEOM in the 100-300 keV Band 



Measured Position 
(Degrees) 



Source Name 



Known Position Position Error 



(Degrees) 



(Degrees) 



(83.59, 22.00) 


CRAB 


(83.63, 22.01) 


0.05 


(201.34,-43.00) 


CENA 


(201.36,-43.01) 


0.03 


(266.25, -30.25) 


IE 1740-29 


(265.96, -29.74) 


0.57 


(268.34, -1.50) 


SWIFT J1753.5-0127 


(268.36,-1.45) 


0.06 


(270.41, -26.50) 


GRS 1758-258 


(270.30, -25.73) 


0.78 


(288.81, 11.00) 


GRS 1915+105 


(288.82, 10.97) 


0.04 


(299.70, 35.00) 


CYGX-I 


(299.59, 35.20) 


0.22 



with variations on monthly timescale s with an average b est fit 
power law spectral index of F = 2.13 (iKataoka et alj|2010l) . The 
LAT data have also been fit using a cut-off power law model re- 
sulting in a better ;|f2 with F = 2.07 and^c = 42.2 GeV. The 
presence of radio synchrotron emission suggests that relativistic 
electrons must be present and must inverse Compton scatter the 
cosmic microwave backgrou nd and infrared ph oton flux, produc- 
ing a hard X-ray y-ray flux {Abdo et al."2009'). Fig. |6la) shows 
the LAT power law spectrum from Kataoka et al. (2010) extrap- 
olated down to GBM energies. The LAT power-law spectrum ex- 
tended down to 12-25 keV falls well below the INTEGRAL data 
points but fits the GBM data (which lie nearly an order of mag- 
nitude below the INTEGRAL points) reasonably well. The GBM 
data do not show any significant variability during the mission so 
far, based on the light curves published on the GBM Occultation 
Project websitqB- 



5.2. 3C 273 

The blight close flat spectrum radio quasar 3C 273 has been 
a frequent target of X-ray and y-ray observations since the 
1970s. Observations with the HEAO A2 experiment in the 2- 



http://heastro.phys.lsu.edu/gbm/ 



60 keV energy rang e are best fit with a power law of F = 1.41 
(IWorrall et al.ll 1979b while HEAO A4 measurement s in the 13- 
180 keV energy range give a slope of F = 1.67 (P rimini et alj 
1979). OSSE observations are well fit by a power law model 
with F - 1.71 and do not show significant spectral variability 
despite flux variabihty d uring observations. With INTEGRAL, 
ICourvoisier eTaP (l2003l) fit the 25-100 keV IBIS/ISGRI spec- 
trum to a power law with F = 1.95 + 0.2 and the 20-200 keV 
SPI spectrum to a power law with F - 1.66 + 0.28. Later 
INTEGRAL and XMM-Newton observations (IChernvakova et all 
2007i) confirmed a softening of the power law spectrum with 
time, reporting a photon index F = 1.82 + 0.01 with observa- 
tions from 2003-2005. Such a soft spectrum can be produced 
by optically thin inverse Compton emission from relativistic 
electrons, but the harder spectrum observed by early experi- 
ments ca nnot be produc ed this way unless, as suggested by 
IChernva kova et all(l2007l) . there is a significant density of pro- 
tons in the y-ray emission region. COMPTEL and EGRET re- 
ported a_brea^_a^~]_MeV with a slope of F - 2.4 above 1 
MeV d Johnson et al.lll995l) . COS B observations at y-ray ener- 
gies (50-800 MeV) show a power law of F = 2.5 (^Bignami et al.l 
1981). However, spectr al hardening was observed at a time when 
the source brightened (Ivon Montignv et al.ll 19971: ICoUmar et al] 
12000 ). Fermi/LAI has observed variations in the photon index 
fromF = 2.4 + 0.2 to 3.3+0.3 and a clear hardening of the spec- 
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trum at times when the flux increases (ISoldiet al.ll2009l) . The 
GBM light curves show distinct variabiHty over the energy range 
12-300 keV. The emission has been modeled by synchrotron 
self Compto n, inverse Compton , and a proton-initiated cascade 
model, and .von Montipnv et al] (Il997h has fit the high-energy 
spectrum satisfactorily with an empirical model of the form 



F(E) = N 



(E/Eb)-^' 
1 + (E/EbY 



(2) 



FiglSb) shows the 4-year average G BM spectrum together with 
the LAT and COMPTEL data and the lvon Montignv et alj ( 1 19971) 
model with Fi = 0.7, r2 = 0.8, and Eb = 2.36 MeV with the 
model and data agreeing up to ~ 10 GeV. 



a timescale of ~ 3 years (IWilson-Hodge et ani201 ih . GBM has 
presented an average spectrum best fit by a broken power law 
model with Fi = 2.06 ± 0.01, F2 ^ 2.36 ± 0.05 and a break 
energy at 98 ± 9 keV (fwilson-Hodge et alJl2012h . INTEGRAL 
(iJourdain & RoquesI |2009|) fit ~ 5.5 years of SPI data to a 
broken power law spectrum covering 20 keV to 6 MeV with 
slopes Fi - 2.07 and F2 - 2.23 and a break energy fixed at 100 
keV. When the break energy is left as a free parameter, Fi = 
2.04 and F2 = 2.18 with a break energy at 62 keV. Fig. gjd) 
shows the INTEGRAL (red dashed) and GBM (solid black) mod- 
els extended to 10 GeV together with results from COMPTEL 
and Fenni/LAiT. The GBM broken power law spectrum provides 
excellent agreement up to 100 GeV. 



5.3. Cen A 

Cen A is one of the br ightest AGN detected f rom radio energies 
up to the TeV range (! Aharonian et al.' '2009") and in ultra-high 
energy cosmic rays by Auger (JAbrahamet al. 2009). Gamma- 
rays are produced by the jet in the central core and by inverse 
Compton scattering of microwave and infrared-optical photons 
in the giant radio lobes. At hard X-ray energies, Ginga and 
balloon-flight observations revealed a power law fit with a sloge 
of F ~ 1.8 with a possible break at ~ 180 keV (Miyazaki et al] 
Il996 ). OSSE, COMPTEL, and EGRET observations during 
an intermediate emission state showed a three-segment broken 
power law with breaks at Ei = 150 keV and £2 - 16.7 MeV 
and slopes of Fi - 1.74, F2 = 2.3 and, F3 - 3.3, while a fit dur- 
ing a low emission state showed breaks at E[ = 140 keV and 
£2 ^ 590 keV a nd slopes of Fi = 1.73, F2 = 2.0 and, F3 = 
2.6 (ISteinleetal.lll998D . BATSE Earth occultation data were 
fit with a s lope of F = 1.84 an d no evidence of a break out 
to 1 MeV dW heaton et aT 1996). INTEGRAL observations by 
iBeckmann et al. (2011) are well fit to an absorbed cut-off" power- 
law model with F = 1.73 and Ec = 434 keV Th e third EGRET 
catalo g reports a power-law fit with F - 2.58 dHartman et al.l 
1 19991). LAT measu rements show a power-law fit with F - 2.67 
dAbdo et al.l 120101. consistent with previous measurements at 
these energies. The observations from LAT energies through the 
radio can be fit with a synchrotron/s ynchrotron self-Co mpton 
model with a single emission region dAbdo etal.1 1201 Ol) or by 
a combined lepto-hadronic model in which the hard X-ray /soft 
y-ray emission is produced by a combination of electron syn- 
chrotron radiation, proton syn chrotron, and inverse Compton 
emission dRevnoso et alj|201l]) . In Fig. |6tc) the Reynoso et al. 
lepto-hadronic model fits well with the LAT, COMPTEL, and 
GBM data over the range from 12 keV to 10 GeV. 



5.4. Crab 

Although the Crab has long been considered to be a steady 
"standard candle", variability was seen with several instruments 
on board CGRO with COMPTEL and EGRET showing vari- 
ations in the 1 MeV to 150 M eV range (De Jager et al. 1996) 
and BATSE ( Ling & Wheatonl2003h suggested flux and spectral 
variability in the 35-300 keV energy range on weekly timescales, 
but showed long-term average spectra consistent with a broken 
power law with F ~ 2.1 and F ~ 2.4 and a break at ~ 100 keV, 
though the BATSE short-term variability results were not con- 
firmed. GBM, INTEGRAL, RXTE, and Swift data have shown 
that the 15-100 keV flux varies by several percent per year on 



6. Conclusions and Future Work 

All-sky images of the hard X-ray /soft y-ray sky have been gen- 
erated with IDEOM using GBM CTIME data. Th ese are the 
first a ll-sky images using Earth occultation since IShaw et aP 
(12004 . The maps cover the 12-50 keV, 50-100 keV, and 100-300 
keV energy bands and span ~ 1500 days. Roughly 45 known 
sources have been identified by the IDEOM imaging. Also 16 
sources have been added to the GBM EOT catalog through 
cross-correlating features in maps with sources in the ^wZ/f/BAT 
and INTEGRAL catalogs. Most of the calculated source positions 
from IDEOM are within ~ 0.5° of the known source position, 
which is the minimum angular resolution as constrained by the 
Earth's atmosphere. No sources were detected in the 100-300 
keV image that were absent from the GBM-EOT catalog down 
to ~ 32 mCrab. This result suggests that the discrepancy between 
the MSEC and JPL BATSE occultation analyses ( Har mon et al.l 
2004) is not due to unmodeled sources, but rather that a more 
likely cause of the difference is the background model. Future 
work with IDEOM includes updating the current images with 
more recent GBM data to increase sensitivity and continuing to 
search for additional sources, and improving IDEOM's search 
algorithm. 

Spectral analysis of the four persistent sources significantly 
detected by both GBM and LAT (NGC 1275, the Crab, 3C 
273, and Cen A) has been performed using GBM, INTEGRAL, 
COMPTEL (where available), and LAT data. For NGC 1275, 
the extrapolation of the average Fenni/LALT spectrum down to 
GBM energies shows results consistent with GBM-EOT. For 
3C 273 an d Cen A, the GB M data are again consistent with 
model fits (Ivon Montignv et al. 1997; Reynoso et al. 2011) ex- 
te nded down from LAT ener gies. For the Crab, the fit suggested 
in lWilson-Hodge et alj d2012.) has been extrapolated up to LAT 
energies and again provides a good fit. 
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Figure 3. All-sky image in the 12-50 keV band for 2008 August 12 to 2012 September 14. Sources marked with blue triangles 
denote sources added to the GBM catalog based on imaging work. Contours have been over plotted at 3. 5, 5, 10, and 20 sigma. 
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Figure 4. All-sky image in the 50-100 keV band for 2008 August 12 to 2012 September 14. Contours have been over plotted at 3. 

5, 5, 10, and 20 sigma. 
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Figure 5. All-sky image in the 100-300 keV band for 2008 August 12 to 2012 September 14. Contours have been over plotted at 3. 

5, 5, 10, and 20 sigma. 
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Figure 6. INTEGRAL (red squares), GBM (black diamon ds), COMPTEL (gree n asterisks) LAT (magenta triangles) (a) NGC 1275 
spectrum with INTEGRAL, GBM, and, LAT data and the iKataoka et alj ( 120101) power law spectrum extrapolated from LAT ener- 
gies down to 12 keV. (b) 3C 273 spectru m with model of the fo rm presented ' von Montigny et al] al997i) overplotted. (c) Cen A 
spectrum with lepto-hadro nic model from iRevnoso et al.l (1201 ih shown. (d)Crab spectrum w ith GBM 3 year average spectrum 
dWilson-Hodge et alj I20T2I) overplotted in (black) solid line, and INTEGRAL/SPl spectrum dJourdain & Roquesll2009h in (red) 
dashed line. 
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